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SPECIFIC COOLING CAPACITY OF LIQUID NITROGEN 


Robert A. Kilgore and Jerry B. Adcock 
Langley Research Center 

ABSTRACT AND SUMMARY 

When cryogenic nitrogen wind tunnels are cooled by injecting liquid 
nitrogen directly into the tunnel, the specific cooling capacity of the 
nitrogen consists of the heat absorbed in vanning and vaporizing the liquid 
pl\is the heat absorbed in warming the gaseous nitrogen to the tunnel stagna- 
tion temperatiu’e. The specific cooling capacity of nitrogen has been calculated 
for a simplified model based on this method of cooling by tising a National 
Bureau of Stsmdards program for the thennodynamic properties of nitrogen and the 
results fitted with a relatively simple equation having tunnel stagnation pres- 
sure and stagnation temperature as the independent variables- This report 
describes the assumed cooling process, describes the method used to calculate 
the specific cooling capacity of liquid nitrogen, gives the sinqple equation 
fitted to the calculated specific cooling capacity data, and presents in 
graphical form calculated values of the specific cooling capacity of nitrogen 
for stagnation temperatvires from saturation to 350 K and stagnct^on pressvires 
from 1 to 10 atmospheres. 


INTRODUCTION 


Based on theoretical studies and experience in the United States 
with the Langley 0. 3-meter transonic cryogenic tunnel, the cryogenic 



wind tunnel concept, using nitrogen as the test gas, has been shown to offer 
many advantages with respect to attaining full-scale R^nolds number at reason- 
able levels of dynamic pressure in tunnels of moderate size. * Because of 
the advantages offered by cryogenic oi>eration, studies are underway for two 
large fan-driven cryogenic tiannels. 

A C3ry genic transonic tunnel, to be known as the National Transonic 
Facility (NTF) and to be located at the NASA Langley Research Center , will 
provide test Reynolds numbers up to 120 million* which is an order of magni- 
tude increase in Reynolds number capability over existing U.S. tunnels. 
Preliminary studies are also being made for a cryogenic transonic tunnel 
to meet the requirements of the European NATO countries. The proposed 
Large Eurox>ean High Reynolds Number Tunnel (LEHRT) will provide test Reynolds 
numbers of at least UO million at transonic speeds.^ 

In addition to the usual aerodynamic and structiural aspects of the 
wind tunnel, the designer of a cryogenic wind tunnel must also provide for 
the cryogenic systems which are unique to this new tunnel concept. In order 
to design properly the cryogenic systems associated with a cryogenic wind 
tunnel, it is necessary to 1l.iow fairly accurately the amoimt of liquid 
nitrogen that will be required for both cool-down and operation of the 
tunnel. Also, since a major part of the direct operating cost of a cryogenic 
tunnel is the cost of the liquid nitrogen, it is essential to know the 
liquid nitrogen requirements in order to properly assess the operating costs. 


^Reynolds numbers quoted are based on a reference length of 0.1 times 
the square root of the test section area. 
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Experience with the Langley C. 3-meter transonic cryogenic tunnel has 
shown that the amount of liquid nitrogen required for cool-down of the 
tunnel structure is predicted adequately by the method of Jacobs. 

The heat to be removed from the tiumel sifter cool-down and while it 
is running consists of the heat conducted through the walls of the tunnel 
from the surroundings, internal sources of heat such as heated actuators 
or internal strain gage balances, and the heat added to the tunnel circuit 
by the drive fain. With a knowledge of the design of the tunnel and the 
properties of the materials used in its construction, the heat conducted 
through the walls can be calculated with sufficient accuracy in a straight- 
forward manner using standard heat transfer equations. Great acciiracy in 
calculating the heat conducted through the walls is generally not required 
since, typically, the heat conducted through the walls amounts to only 
about 10 percent of the total heat to be removed at subsonic speeds and 
about 1 percent at trsmsonic speeds.^ The heat to be removed due to 
internal sources such as heated actuators or balances is generally small 

O 

enou^ to be negligible. If, under special circumstances, heat from these 
intenial sources must be known, it may be determined experimentally or 
estimated from a detailed emalysis of each source. 

Calculating the heat added by the drive fsui can also be done in a 

9 

straightforward manner by using the method of Adcock and Ogbxurn who, by 
taking into account the real gas properties of nitrogen, have determined 
the drive power required for Isentropic compressions of cryogenic nitrogen. 
Their analysis, which covers stagnation temperatures from 80 to 310 K and 
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stagnation pressures from 1 to 8*8 atmospheres, shows the power required 
to cosq)re8S cryogenic nitrogen to be less than that calculated for an 
ideal diatomic gas by as much as 9*^ percent at the highest pressure* This 
reduction is, of course, in addition to the basic reduction in drive power 
with decreasing temperature, that is, power varies as the square root of 
the stagnation temperature.^ The reader is thus advised to consult refer- 
ence 7 for the proper method of calcxilating the heat added to the stream 
by the drive fan. 

Once the total heat load is calculated, it is necessary to know the 
specific cooling capacity of liquid nitrogen in order to calculate the 
amount of liquid nitrogen required for tunnel operation* Very precise 
values of the specific cooling capacity can be obtained by suitable 
manip\ilation of the data contained in a National Bureau of Standards (NBS) 
computer program based on the work of Jacobsen *^^ However, because of th 
large size of the NBS program it is impractical to use for direct calcula- 
tion of specific cooling capacity with relatively small desk-top program- 
mable computers* Therefore a relatively simple equation which can easily 
be used with such digital computers has been fitted to the specific cooling 
capacity data determined from the NBS program. 

This report describes the assumed cooling process, the method used to 
csilculate the specific cooling capacity of liquid nitrogen from the NBS 
program, and gives a relatively simple equation for calculating specific 
cooling capacity as a function of tunnel stagnation values of temperature 
and pressure. In addition, the specific cooling capacity of nitrogen is 
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given in graphiced form for stagnation temperatures from saturation to 350 K 
and stagnation pressures from 1 to 10 atmospheres. 

SYMBOLS 

h specific enthadpy* kJoules/kgm 

p absolute pressiu'e, atm (l atm * 101.3 TsS/v?) 

q heat, kJoules 

s specific entropy, kJoules/kgm - °K 


T absolute tonperature, K 

B specific cooling capacity of nitrogen, kJoules/kgm 

Subscripts : 

a supercooled liquid (obtained .by isentropic compression 

from p = 1 atm) , p = p.^. 

b liquid at the saturation bovindary, P = p.|. 

b" gas at the satioration boundary, p “ p.j. 


t gas at the tunnel stagna m conditions , p = Pt » T = T.^. 
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SPECIFIC COOLING CAPACITY OF NITROGEN 


Assumed Cooling Process 

Cooling in the 0.3-meter transonic cryogenic tunnel is accomplished by 
injecting liquid nitrogen (I1N2) directly into the tunnel circviit. The LNg 
is stored in a large insulated tank kept at atmospheric pressure. A cen- 
trifugal p\2nq> is used to increase the pressure of the LNg to a pressure 
slightly greater than the tunnel stagnation pressure p^. The LN2 is then 
injected into the tunnel at the diffuser downstream of the test section where 
it vaporizes. The resulting nitrogen gas absorbs additional heat so that 
its final temperature and pressure are identical to the stagnation 
values p^ and T^. Since this simple and direct method of cooling the 
0.3-meter transonic cryogenic tunnel has proven to be satisfactory, the 
specific cooling capacity of nitrogen will be determined based on this 
method. The basic assumptions are thei'efore; 

1) Liquid stored at atmospheric pressure 

2 ) Liquid compressed isentropically to p = 

3) Liquid injected directly into tunnel circuit 

4) Warming and vaporization of liquid occTzrs at p = 

5) Warming of gas to occurs at p = p^ 

The basic assumptions Virch have been made are an oversimplification of the 
actual cooling process. For example, even for unpressxirized LNg storage 
tanks, the head of liquid nitrogen makes the simplifying assumption of 
storage at exactly atmospheric pressure invalid. Also, in practice the 
LNg may at times be pressurized by the centr*' f’-igal pump to a pressure 
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several atmospheres above p^.. The a8S\iB^tlon that warming suad vaporization 
of the liquid nitrogen takes place at a constant pressiure equal to the stag- 
nation pressure, p^, is, of course, a great simplification since in practice 
the liquid is injected in the low-pressure region of the diffuser and the 
actual warming and vaporization process occiirs under conditions of continually 
changing pressure. The simplifying assumptions not only allow the cooling 
capacity to be calculated in a relatively siit^le and straightforward manner, 
but, more importantly, as csm be found from a more rigorous analysis of the 
cooling process, the simplifying assumptions introduce only negligible errors 
in the calculated values. From basic thermodynamic considerations it is 
known that areas in a temperature -entropy diagram represent heat. In order 
to show the various steps of the cooling process, the temperature-entropy 
diagram for the heat absorbing process being considered is shown schematically 
in figure 1. For the process proceeding from point a (supercooled liquid 
P ® Pt *1^® "to isentropic compression from p = 1 atm) to point t (gas 
at p^ and T^ ) , the heat absorbed by a unit mass cf nitrogen is propor- 
tional to the area bounded by the curve a b b' t, the entropy axis, and 
the vertical lines at s^ and s^. The total heat absorbed per unit mass, 
defined as the specific cooling capacity 3, is proportional to the sum of 
the three areas 

® %-b' * %'t 

where * heat absorbed in warming the liquid nitrogen 
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= heat of vaporization 


absorbed in warming the gaseous nitrogen 

Cooling Due to Heat Absorbed in Warming the 
Supercooled Liquid Nitrogen 

The area identified as q^^^ in figure 1 is proportional to the heat 
absorbed in warming a unit mass of the liquid nitrogen from approximately 
the storage temperature at atmospheric pressure to the saturation temperature 
corresponding to the stagnation pressure p^. By assuming that the liquid 
nitrogen is compressed isentropically from its storage pressure to a pres- 
sure equal to p^, the calculation of cooling capacity is greatly simplified* 
’Phis results from the fact that the heat absorbed between any two points 
lying on the same isobar is equal to the difference in the specific enthalpy 
between the two points, [it should be noted that the Isobars are so com- 
pressed in the liquid region at the saturation bound€u*y that the temperature 
rise due to t'le isentropic compression is very smeill) being only 0.2U^K for 
pressurization from 1 atm to 10 atm.) 

Thus, the heat absorbed in wanning the supercooled liquid can be 
calculated as 


'lab " ^b ~ \ 

where h = specific enthalpy 


a = supercooled liquid, p = 


b = liquid at the saturation boundary, p * 
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Since the hoiling point of liquid nitrogen Increases with increasing 
pressure, the advantage of storing the liquid nitrogen at atmospheric 
pressure rather than at elevated pressure is apparent from figure 1 vhere 
it can he seen that the heat absorbing capacity due to warming a unit mass 
of the liquid, represented by the area decreases as the long-term 

storage pressure, and hence the equilibrium storage temperature,, is 
increased. The decrease in cooling capacity with long-term storage 
pressure will be described in more detail in a later section. 

Cooling Due to Vaporization of the Liquid Nitrogen 

The aiea q^^j- in figure 1 is proportional to the heat absorbed per 
unit mass of nitrogen in the change of phase from liquid at b to gas 
at b**. Again, since this process is assumed to take place along the pt 
isobar, csin be caJ.culated directly from the difference in the 

specific enthali^ of the two phases across the saturation boundary. 

%b' “ \ 

where b = liquid at the saturation bo\mdary, p = p,j. 

b' = gas at the satviration boundary, p ■ 

Cooling Due to Keat Absorbed in Warming the Gaseous Nitrogen 

The area in figure 1 is proportional to the heat absorbed as 

a unit mass of nitrogen gas is warned from the saturation temperature 
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corresponding to to the tunnel stagnation temperatixre . Since this 

wELnning process is assumed to take place along an isobar (p Pt)* 
heat absorbed per unit mass is again equal to the d-fferences between the 
specific enthalpies of the gas at the end points. Thus 




(3) 


where = gas at the saturation boundary, p * 


t = gas at the tunnel stagnation conditions, p - p^> T = 


Total. Specific Cooling Capacity 


The three areas in the temperature-entropy diagram of figure 1 were 
considered separately in the preceding sections in order to describe the 
three separate heat absorbing processes contributing to the total specific 
cooling capacity. However, because the entire cooling process is assumed, 
for simplicity, to take place along the p.^ isobar, the total heat 
absorbed depends only on the difference between the specific enthalpy of 
the nitrogen at the end points of the heat absorbing process, a ar.d t. 
That this is true can also be seen by adding the sepo 2 *ate contributions 
from equations (l), (2), and (3) to obtain the total specific cooling 
capacity 


^ = ‘lab %h' * %'t “ K - - l»b) + K " 

or 

^ ^ \ p = Pt 

^Failure to properly evaluate q, resulted in the values of the 
specific cooling capacity for nitrogen given in reference T being too large 
by about 10 percent at the highest temperatures. 
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BQUATIOH FOR SPECIFIC COOLING CAPACITY AS A FUNCTION 
OF AND T|. 

Values of the specific booling capacity of nitrogen, 3* were calcu- 
lated from the National Bureau of Stemdards counter program for pressures 
from 1 tolOatm and for temQreratures from saturation to 350 K. The 
calculated values of 3 were then fitted witn an equation of the form 

e ■ ®o ’’t ♦ »i ♦ % I- 

^t 

Where Bq = 1.0379 - 3.9157 x 10"^ p^ 

= 1.2125 X 10^ + 2.1577 Pt 
Bg =* £^.6585 X 10^ - 3.9122 x lO^ p^ 

T^ = stagnation temverature, K 
p. stagnation pressure, atm 

u 

The dimensions of 3 calculated from equation (5) are kJ/kgm. The values 
of 3 calculated from equation (5) agree well with the vedues of specific 
cooling capacity o^ nitrogen calculated using the NBS program and the 
method describe i In the preceding section (eq. (4)) for absolute pressures 
from 1 to 10 atm and absolute temperatures ^rom saturatior ‘-o 350 The 
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arera^ error in 0 is less tbsn 0.3 percent vith tbe iwurimHn error telog 
0.9 percent for values of 0 calculated at the Mgher pressures vith 
at the saturation boundary. 

The specific cooling edacity of liqt^d nitrog«i« 0* as a function of 
tunnel stagnation teqterature for various val'tes of tunnel stagnation pres>' 
sure uas calculated using equation (3) and is presented graphically in 
figure 2 The conditions covered are likely to cov^ all practical wind 
tunnel designs* naaely, stagnation teaperatiires froa saturation to 350 K 
and stagnatiem jiressures froa 1 to 30 ata. As can be seen froa figure 2* 
the specific cooling capacity of liquid nitrogen is a weak function of 
stagnation pressure and a fairly strong nearly linear function of stagna- 
tion teoperature. 

As previously noted, the cooling capacity of liquid nitrogen is 
reduced if the liquid is stored at elevated pressiires and the bulk 
teoperature of the liquid rises. The decretise in cooling capacity due 
to storage of the liquid nitrogen at elevated pressiires is shown in fig- 
ure 3 as a function of storage pressure. The values of A0 are based 
on the assimption that the t«q>erature of the liquid is eqiml to the 
saturation toiperature corresponding to the storage pressure. For long- 
term storeige at pressvires above 1 atm the values of 0 given by equa- 
tion ( 5 ) or in figure 2 are reduced by the value of A0 given in 
figure 3 for the actual storage pressure. 

As an exEunple of the reduction in cooling capacity due to storage 
at elevated pressures we can consider cooling a tunnel with stagnation 


12 



conditions of • 5 atm and ■ 100 K. Fron figure 2 (or eq. (5)) ve 
find 6 • 215.2 kJ/k^ based on storage at 1 atm (T > 77>35 K) vith 
isentropie compression to 3 atm (resulting in a siqpercooled liqjuld at 
T • 77.4$ K). Had the liquid nitrogen been stored at 7 atm ve find 
from figure 3 that the reduction in cooling edacity, is 4$.0 
Under thr conditions, cooling capacity is reduced by 21$. 

CQBCUJDIHG RBIABKS 

The specific cooling capacity, 6, of nitrogen has been calculated 
for a simplified cooling model based on the cooling method developed for 
the Langley 0.3^ transonic cryogenic tunnel. Calculations of $ were 
made by using a HationcLl ftxreau of Standards program for the themK>dynamic 
properties of nitrogen and the results fitted with a relatively sisqale 
equation of eight coefficients having tunnel stagnation pressure and 
stagnation tanperature as the independent variables.* The eq\iation can 
be used with relatively small desk-top coii^uters to calculate the specific 
cooling capacity of nitrogen for tiumel stagnation temperatures from 
saturation to 330 K and tunnel stagnatic.i pressures from 1 to 10 atmospheres. 
The values of 3 calculated from the relatively siople equation generally 
agree with the more precise calciilations made directly from the National 
Bureau of Standards program to within 0.3 percent with a maximum differ- 
ence of only 0.9 percent for values c€d.culated with tunnel stagnation 
temperature neeu* the saturation boundary. 
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Figure L* A schematic temperatura-entropy diagram for the heat 
absorbing process being considered. 
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Storage Pressure, ATM 

Figure 3.- Reduction in cooling capacity of liquid nitrogen due to storage 
at pressures greater than I atm. 















